A Robust GPS Navigation Filter Based on Maximum Correntropy Criterion with Adaptive Kernel Bandwidth

Multiple forms of interference and noise that impact the receiver’s capacity to receive and interpret satellite signals, and consequently the preciseness of positioning and navigation, may be present during the processing of Global Positioning System (GPS) navigation. The non-Gaussian noise predominates in the signal owing to the fluctuating character of both natural and artificial electromagnetic interference, and the algorithm based on the minimum mean-square error (MMSE) criterion performs well when assuming Gaussian noise, but drops when assuming non-Gaussian noise. The maximum correntropy criteria (MCC) adaptive filtering technique efficiently reduces pulse noise and has adequate performance in heavy-tailed noise, which addresses the issue of filter performance caused by the presence of non-Gaussian or heavy-tailed unusual noise values in the localizing measurement noise. The adaptive kernel bandwidth (AKB) technique employed in this paper applies the calculated adaptive variables to generate the kernel function matrix, in which the adaptive factor can modify the size of the kernel width across a reasonably appropriate spectrum, substituting the fixed kernel width for the conventional MCC to enhance the performance. The conventional maximum correntropy criterion-based extended Kalman filter (MCCEKF) algorithm’s performance is significantly impacted by the value of the kernel width, and there are certain predetermined conditions in the selection based on experience. The MCCEKF with a fixed adaptive kernel bandwidth (MCCEKF-AKB) has several advantages due to its novel concept and computational simplicity, and gives a qualitative solution for the study of random structures for generalized noise. Additionally, it can effectively achieve the robust state estimation of outliers with anomalous values while guaranteeing the accuracy of the filtering.


Introduction
Due to multipath interference and non-line-of-sight (NLOS) reception, global navigation satellite systems (GNSS), such as the Global Positioning System (GPS), generate significant errors in urban settings [1,2].Due to diffraction from adjacent structures like the ground and water surfaces, structures, moving cars, hills, trees, etc., and multipath effects arise when GPS signals travel via many pathways before reaching a receiver site.Numerous estimate strategies have been investigated to remove the location inaccuracy brought on by multipath.Although multipath errors are a class of uncorrelated mistakes not eliminated by observation differentiation, they primarily restrict the performance of high-precision GPS receivers.The strategies for interference minimization are among the most crucial aspects of enhancing the efficiency of the GPS system.
The outliers are assessments that considerably differ from the majority of the observation sample.The irrelevancy of distinct instances and the emergence of uncertainty advantages due to its novel concept and computational simplicity.Also verified from the simulation results is a novel generic solution for the study of random structures with the generalized noise.In order to enhance the performance of EKF applications when non-Gaussian disturbances exist, MCCEKF uses the MCC rather than the MMSE as the optimization criterion for GNSS navigation processing.With the observed distance and its variance, regulated to eliminate the effects of colored Gaussian noise, the range error space may be restricted.The suggested novel MCCEKF-AKB has been shown to be more efficient at handling non-Gaussian data by adding an adaptable kernel bandwidth component to the MCCEKF.
The remainder of this paper is organized as follows.A brief review of the maximum correntropy criterion is reviewed in Section 2. In Section 3, the basic description of the MCC-based EKF (MCCEKF) is presented.In Section 4, the MCCEKF with adaptive kernel bandwidth (MCCEKF-AKB) is introduced.The proposed MCCEKF-AKB's performance compared to the EKF and MCCEKF techniques are assessed using illustrative examples based on simulation experiments in Section 5. Finally, conclusions are given in Section 6.

Maximum Correntropy Criterion
Considering the autocorrelation of random phenomena, correntropy was first introduced to quantify consistency across delays.Later, it was expanded to evaluate the localized similarity of any two random variables.The correntropy between X and Y can be defined as follows: where E[•] represents the expectation function, the kernel function, and the joint probability density function (pdf) of X and Y are denoted as κ σ (X, Y), and f XY (x, y), respectively.In general, the kernel size in terms of the Gaussian function can be expressed as Generally, f XY (x, y) is calculated with a limited availability of data, which leads to computational difficulty.Therefore, the mean of a sample is used to compute the correntropy.
where (x i , y i ) N i=1 are the N samples drawn from f XY (x, y).After taking the Taylor series expansion of the Gaussian kernel function, we have In order to represent the high-order information within the data, the correntropy is a weighted sum for all the even-order moments of the variable (X − Y).Robust state es- timation can be accomplished using correntropy as an optimization strategy because it measures if the degree of similarity exists between X and Y.
Although, several definitions were there, Renyi and Shannon have been created for diverse purposes since Shannon first proposed entropy in 1948 [29].Correntropy, for instance, is frequently used in adaptive filtering and machine learning to determine the local similarity of random variables.The correntropy of two random variables that are presented in the kernel area is defined by information-theoretic and kernel techniques as a generalized local matching score.Several applications have been accepted by the MCC.For instance, machine learning, pattern recognition, and signal processing have all helped the MCC.The MCC and extended Kalman filters are introduced in this section to cope with non-white Gaussian noise, particularly the outliers.

Extended Kalman Filter with Maximum Correntropy Criterion
Considering a model equation in the form of one discrete-time is where x k ∈ n and w k ∈ m denote the nonlinear state vector and the process noise vectors, respectively, and z k ∈ m and v k ∈ m represent the measurement vector and the measurement noise vector, respectively.Q k is the process noise covariance matrix where R k is the measurement noise covariance matrix.
The vectors w k and v k in Equations ( 1) and ( 2) are zero-mean Gaussian white sequences having zero cross-correlation with each other: where the superscript "T" denotes the matrix transpose.The symbol δ ik stands for the Kronecker delta function given by The details algorithm for discrete-time EKF are as follows: (1) Initialization of state vector x0 and state covariance matrix P 0 .
(2) Predictions of the state vector and state covariance matrix are (4) Updating of the state vector: (5) Updating of the error covariance: where I is the unity matrix.For a linear dynamic system represented by (for Equation ( 5)), the prediction of the state vector is then written as xk|k−1 = Φ k−1 xk−1 (for EKF algorithm step 2) and no Jacobian calculation is required.The relations between the linear approximation equation and the measurement matrices of the system can be expressed as To obtain the filter gain K k for a Kalman filter, the performance index is minimized in terms of the MSE condition: where tr(•) denotes the trace of a matrix.In the absence of uncertainty in the process and measurement noise covariances, the performance index J MSE attains the global minimum by using the standard Kalman filter, taking the partial derivative of P k with respect to K k , using where δx k is the state prediction error as δx k = xk|k−1 − x k and we have The linear Kalman filter's nonlinear variant can be obtained in detail from the linear model in [9] using the MCC as the base.Both of the two cost functions deal with various statistical parameters of the random variables.The cost function based on the MSE criterion is given by Equation (7), which is a function of the covariance matrix of the estimation error.On the other hand, the cost function based on the MCC given by Equation ( 11) is a function where the quantity involved is the innovation sequence z k − h(x k|k−1 ) weighted by R −1 k and the residual xk − xk|k−1 weighted by P −1 k|k−1 .The weighted least squares (WLS) method can be used to create an alternate objective function and improve its robustness.Consequently, the objective purpose of MCC might be described as where ||•|| A represents the weighted norm of A. Taking the derivative of the objective function J MCC with respect to xk and setting its value to zero, and the following equation can be obtained: which can be written as where and Sensors 2023, 23, 9386 6 of 17 We can write then, finally, the factor can be obtained.
The MCC based EKF results for the estimation with the updated gain can be written as which can also be expressed as Figure 1 provides the flow chart for the maximum correntropy criterion-based extended Kalman filter (MCCEKF).The computation procedures are summarized as follows: (1) Initialization of the state vector and state covariance matrix: x0 and P 0 ; (2) Prediction of the state vector and the state covariance matrix: xk|k−1 and P k|k−1 ; (

MCCEKF with Adaptive Kernel Bandwidth
Considering the initial measurement k z and the optimal prediction

MCCEKF with Adaptive Kernel Bandwidth
Considering the initial measurement z k and the optimal prediction xk|k−1 acquired in the previous step, the innovations sequence is given by The innovation shows the difference between the predicted measurement ẑk|k−1 and the actual measurement z k , while z k gives the extra data that the filter has access to as a new observation of the result.The weighted innovation, K k (z k − ẑk|k−1 ), acts as a correction to the predicted estimation xk|k−1 with the estimated xk .

Innovation Information for Failure Detection and Adaptive Algorithms
Further, fetching both side variances, the obtained theoretical covariance, and the covariance matrix of the innovated sequence can be expressed as which can also be rewritten as Given that the length of the window is one and the degree of divergence (DOD) [15] is used as the record of the innovation covariance matrix, This variable can be adjusted for the adaptive filtering or utilized to identify the divergence or outliers.Alternately, the concepts behind Equations ( 26) and ( 28) can be used to establish the parameters for determining the level of change in vehicle dynamics.The other DOD measure, which is frequently employed as a straightforward test statistic for failure identification, is based on the ratio Similar to the above-mentioned parameters, the normalized innovation squared utilized in the case of fixed kernel bandwidth can be given by (υ , where a factor β can be incorporated to obtain the extended version

MCCEKF Based on Adaptive Kernel Bandwidth Mechanism
Employing these upper limits as an outline, the size of the kernel must be adjusted.Additionally, the kernel size is modified in the following way in the measurement-specific processing for outliers: where λ j,k denotes the adaptive factor for the jth measurement element at time step k, z j,k , with the maximum kernel size σ max = 150.
Further, λ j,k must be suitably adjusted for varying noise intensities.The correlation between the innovation factor and its covariance matrix C υ k is the first point to consider in finding the presence of measurement outliers, and α j,k can be defined as where C υ j,k is the jth diagonal element of the innovation covariance matrix with υ j,k being the innovation term υ k for the jth element, and υ j,k = z j,k − ẑj,k|k−1 .The λ j,k is obtained by using the relation Again, Now, we have to find σ j,k based on the adaptive factor λ j,k according to α j,k .
The kernel size can be obtaining once the adaptive factor λ j,k , is presented in Figure 2. Figure 3 shows the basic flow diagram to obtain the prefect kernel bandwidth.Further, the logic for the kernel bandwidth adapting algorithm is described in case studies.
Case 1: If <α j,k > 1>, then <λ j,k is increased> (this indicates that if C υ k is larger than υ 2 j,k , then the kernel bandwidth increases for maintaining the optimal performance).Case 2: If <α j,k < 1>, then <λ j,k is decreased> (this indicates that if C υ k is less than υ 2 j,k , then the kernel bandwidth decreases for maintaining the robustness performance).

( )
, , 1 exp Flow chart for calculation of the kernel bandwidth.
The flow chart for the MCCEKF-AKB is presented in Figure 4, and the algorithm steps are summarized below:  The flow chart for the MCCEKF-AKB is presented in Figure 4, and the algorithm steps are summarized below: (1) Initialization of the state vector and state covariance matrix: x and P ; The flow chart for the MCCEKF-AKB is presented in Figure 4, and the algorithm steps are summarized below: (1) Initialization of the state vector and state covariance matrix: x0 and P 0 ; (2) Prediction of the state vector and the state covariance matrix: xk|k−1 and P k|k−1 ; (3) Obtaining σ j,k based on the adaptive factor λ j,k according to α j,k ; , ma

Results and Discussion
In order to demonstrate the efficiency of the suggested methodology for the context of a time-varying environment, computational simulations were carried out to assess the performance of the recommended strategy to traditional methods of GPS navigation processing.The authors designed the various navigation filters, such as the EKF, MCCEKF, and MCCEKF-AKB for the nonlinear systems, using MATLAB 2022 ® codes and the commercially available programs Inertial Navigation System (INS) Toolbox [30] and Satellite Navigation (SatNav) Toolbox [31], both provided by GPSoft LLC, Athens, Ohio, USA.The proposed filtering algorithm does not affect the navigation signal since there is no information (e.g., velocity) aiding feedback to the tracking loop.The INS Toolbox is used for generating the vehicle's trajectory, while the SatNav Toolbox is used for generating the satellite orbits and pseudoranges.
The satellite skyplot for the seven GPS satellites (the red colored dot followed by the GPS ID number), and the test trajectory of the simulation are shown in Figure 5a,b, respectively.All the satellites are reachable, which is advantageous for ensuring proper navigation satellite geometry and accurate positioning results.The vehicle taken for simula-

Results and Discussion
In order to demonstrate the efficiency of the suggested methodology for the context of a time-varying environment, computational simulations were carried out to assess the performance of the recommended strategy to traditional methods of GPS navigation processing.The authors designed the various navigation filters, such as the EKF, MCCEKF, and MCCEKF-AKB for the nonlinear systems, using MATLAB 2022 ® codes and the commercially available programs Inertial Navigation System (INS) Toolbox [30] and Satellite Navigation (SatNav) Toolbox [31], both provided by GPSoft LLC, Athens, Ohio, USA.The proposed filtering algorithm does not affect the navigation signal since there is no information (e.g., velocity) aiding feedback to the tracking loop.The INS Toolbox is used for generating the vehicle's trajectory, while the SatNav Toolbox is used for generating the satellite orbits and pseudoranges.The satellite skyplot for the seven GPS satellites (the red colored dot followed by the GPS ID number), and the test trajectory of the simulation are shown in Figure 5a,b, respectively.All the satellites are reachable, which is advantageous for ensuring proper navigation satellite geometry and accurate positioning results.The vehicle taken for simulation is assumed starting from the at 121.7775 degrees (east longitude) and 25.1492 degrees (north latitude) and is placed at a location with an altitude of 100 m, which is at −3042329.20 4911080.20 2694074.30 T m at the WGS-84 ECEF coordinate frame.The vehicle is the starting location and is considered as the (0,0,0) point in the local tangent East-North-Up (ENU) frame.
i -th satellite's position in three dimensions, c is the speed of light and b t is the re- ceiver clock offset from system time, and i v  is the pseudorange measurement noise.In order to evaluate the performance when dealing with outlier-type multipath interferences, several additional and randomly produced errors are purposely introduced into the GPS pseudorange observed data while the vehicle is in motion.Since the errors caused by pseudorange observations indicate the accuracy of the measurement designs in the filters, we can concentrate on investigating these errors.The receiver thermal noise, multipath, ionospheric delay, and tropospheric delay are the error sources that affect the GPS outputs.The majority of inaccuracies are subsequently addressed by using the differential GPS (DGPS) mode, although the effects of multipath and the receiver's thermal noise cannot completely be avoided in the current investigation.The error brought on by the thermal noise in the receiver appears as a Gaussian sequence with a mean of 0 m and a variance of 1 m 2 .For instance, using GNSS data to navigate in urban regions or environments with NLOS and multipath receptions is challenging.Additionally, the GPS sensor's ability to locate itself precisely is hampered by the addition of incorrect data caused by these incidents.The additional interference sequences due to the outliers for the pseudorange observations throughout the process are shown in Figure 6.There are 30 outliers in this simulation.During the coding process, these measurement errors could range from a few tens to several hundred meters.The dynamic process of the GPS receiver is a linear one represented by the PV (Position-Velocity) model [3], where three position states, three velocity states, and two GPS receiver clock states (clock bias and drift) are involved, so that the state to be estimated is an 8 × 1 vector.The measurement model is a nonlinear system model, involved in the GPS pseudorange observables for n satellites where (x, y, z) denotes the user position in three dimensions, (x i , y i , z i ) denotes the i-th satellite's position in three dimensions, c is the speed of light and t b is the receiver clock offset from system time, and v ρ i is the pseudorange measurement noise.
In order to evaluate the performance when dealing with outlier-type multipath interferences, several additional and randomly produced errors are purposely introduced into the GPS pseudorange observed data while the vehicle is in motion.Since the errors caused by pseudorange observations indicate the accuracy of the measurement designs in the filters, we can concentrate on investigating these errors.The receiver thermal noise, multipath, ionospheric delay, and tropospheric delay are the error sources that affect the GPS outputs.The majority of inaccuracies are subsequently addressed by using the differential GPS (DGPS) mode, although the effects of multipath and the receiver's thermal noise cannot completely be avoided in the current investigation.The error brought on by the thermal noise in the receiver appears as a Gaussian sequence with a mean of 0 m and a variance of 1 m 2 .For instance, using GNSS data to navigate in urban regions or environments with NLOS and multipath receptions is challenging.Additionally, the GPS sensor's ability to locate itself precisely is hampered by the addition of incorrect data caused by these incidents.
The additional interference sequences due to the outliers for the pseudorange observations throughout the process are shown in Figure 6.There are 30 outliers in this simulation.During the coding process, these measurement errors could range from a few tens to several hundred meters.

Scenario 1: Pseudorange Observable Errors Based on Gaussian Mixture Distribution
Scenario 1 describes how to reduce outlier-type multipath problems that affect the pseudorange observable.Further randomly produced errors have been deliberately incorporated into the GPS pseudorange observational information over an aggregate of five time periods while the vehicle is traveling.
The majority of the receiving-dependent errors are assumed to be avoided in this work using the differential GPS method, but receiver-independent problems like multipath and receiver measurement noise from thermal sources are not.We believe that the additional errors, such as multipath error and the thermal noise, have the Gaussian mixed distribution type, written in terms of the probability density function as follows: which is equivalent with where  7. Further, in Figure 8, the positioning errors of the ENU frame for EKF and MCCEKF are shown.It can be noted that the convergence behavior of MCCEKF is significantly influenced by the kernel bandwidth, which can be fixed by the Banach Fixed-Point Theorem [7] in order to obtain a fixed-point convergence range.Furthermore, Figures 9 and 10 show the different results for MCCEKF versus MCCEKF-AKB, and for the three approaches (EKF, MCCEKF, and MCCEKF-AKB), respectively.The results convey that the rate of convergence speed increases with adaptive kernel bandwidth.It will take a long time for the algorithm to diverge or converge if the kernel bandwidth is too tiny.

Scenario 1: Pseudorange Observable Errors Based on Gaussian Mixture Distribution
Scenario 1 describes how to reduce outlier-type multipath problems that affect the pseudorange observable.Further randomly produced errors have been deliberately incorporated into the GPS pseudorange observational information over an aggregate of five time periods while the vehicle is traveling.
The majority of the receiving-dependent errors are assumed to be avoided in this work using the differential GPS method, but receiver-independent problems like multipath and receiver measurement noise from thermal sources are not.We believe that the additional errors, such as multipath error and the thermal noise, have the Gaussian mixed distribution type, written in terms of the probability density function as follows: which is equivalent with where σ 1 and σ 2 represent the standard deviations of the individual Gaussian distribution, and the error model contamination perturbing parameter is ε.To illustrate the viability of the suggested algorithm's robustness in the form of non-Gaussian distribution, the concerned requirements are ε = 1/30 (i.e., with approximate 3.33% contamination), and σ 2 = 100σ 1 = 100 × 1, which have been used to simulate the outlier type of multipath errors.The results of positioning errors for different techniques have been verified from Figures 7-10.The position errors corresponding to different values of extended kernel bandwidth factor β with respect to the time (s) are presented in Figure 7. Further, in Figure 8, the positioning errors of the ENU frame for EKF and MCCEKF are shown.It can be noted that the convergence behavior of MCCEKF is significantly influenced by the kernel bandwidth, which can be fixed by the Banach Fixed-Point Theorem [7] in order to obtain a fixed-point convergence range.Furthermore, Figures 9 and 10 show the different results for MCCEKF versus MCCEKF-AKB, and for the three approaches (EKF, MCCEKF, and MCCEKF-AKB), respectively.The results convey that the rate of convergence speed increases with adaptive kernel bandwidth.It will take a long time for the algorithm to diverge or converge if the kernel bandwidth is too tiny.
Faster convergence is guaranteed by a wider kernel bandwidth; however, this frequently results in a low performance in impulsive noise environments.The bandwidth can be manually adjusted or improved using trial-and-error techniques in real-world situations and it is not realistic in real-time implementations.Figure 11 shows the Monte Carlo trials for the kernel size against the time for considering the AKB for the MCCEKF, because the adaptive kernel size is related to the selection of σ max , it has resistance towards the outlier interferences [14].Faster convergence is guaranteed by a wider kernel bandwidth; however, this frequently results in a low performance in impulsive noise environments.The bandwidth can be manually adjusted or improved using trial-and-error techniques in real-world situations and it is not realistic in real-time implementations.Figure 11 shows the Monte Carlo trials for the kernel size against the time for considering the AKB for the MCCEKF, because the adaptive kernel size is related to the selection of max  , it has resistance to- wards the outlier interferences [14].

Scenario 2: Pseudorange Observable Involving Outlier Type of Multipath Interferences with Time-Varying Variance in Measurement Noise
For Scenario 2, we considered the pseudorange observable involving outlier type of multipath interferences with time-varying variance in the measurement noise.Figure 12 represents the information sequence affected due to the outlier interference, which indicates the variation in outliers with changing kernel size and time.For a more quantitative analysis, the results for various positioning errors for the ENU frames are represented in Figures 13-16.The position errors corresponding to different values of extended kernel bandwidth factor  with respect to the time (s) are presented in Figure 13.Furthermore, Figure 14 represents the results for EKF and MCCEKF and Figure 15 shows the results for MCCEKF and MCCEKF-AKB, respectively, with a fixed value of kernel bandwidth.It is worth noting that the kernel bandwidth plays a vital role for the MCCEKF, which leads to the results of [9,11,28].Additionally, the positioning errors among the three approaches are shown in Figure 16.The proposed algorithm for MCCEKF-AKB (Figure 4) achieves a Faster convergence is guaranteed by a wider kernel bandwidth; however, this frequently results in a low performance in impulsive noise environments.The bandwidth can be manually adjusted or improved using trial-and-error techniques in real-world situations and it is not realistic in real-time implementations.Figure 11 shows the Monte Carlo trials for the kernel size against the time for considering the AKB for the MCCEKF, because the adaptive kernel size is related to the selection of max  , it has resistance to- wards the outlier interferences [14].

Scenario 2: Pseudorange Observable Involving Outlier Type of Multipath Interferences with Time-Varying Variance in Measurement Noise
For Scenario 2, we considered the pseudorange observable involving outlier type of multipath interferences with time-varying variance in the measurement noise.Figure 12 represents the information sequence affected due to the outlier interference, which indicates the variation in outliers with changing kernel size and time.For a more quantitative analysis, the results for various positioning errors for the ENU frames are represented in Figures 13-16.The position errors corresponding to different values of extended kernel bandwidth factor  with respect to the time (s) are presented in Figure 13.Furthermore, Figure 14 represents the results for EKF and MCCEKF and Figure 15 shows the results for MCCEKF and MCCEKF-AKB, respectively, with a fixed value of kernel bandwidth.It is worth noting that the kernel bandwidth plays a vital role for the MCCEKF, which leads to the results of [9,11,28].Additionally, the positioning errors among the three approaches are shown in Figure 16.The proposed algorithm for MCCEKF-AKB (Figure 4) achieves a

Scenario 2: Pseudorange Observable Involving Outlier Type of Multipath Interferences with Time-Varying Variance in Measurement Noise
Scenario 2, we considered the pseudorange observable involving outlier type of interferences with time-varying variance in the measurement noise.Figure 12 represents the information sequence affected due to the outlier interference, which indicates the variation in outliers with changing kernel size and time.For a more quantitative analysis, the results for various positioning errors for the ENU frames are represented in Figures 13-16.The position errors corresponding to different values of extended kernel bandwidth factor β with respect to the time (s) are presented in Figure 13.Furthermore, Figure 14 represents the results for EKF and MCCEKF and Figure 15 shows the results for MCCEKF and MCCEKF-AKB, respectively, with a fixed value of kernel bandwidth.It is worth noting that the kernel bandwidth plays a vital role for the MCCEKF, which leads to the results of [9,11,28].Additionally, the positioning errors among the three approaches are shown in Figure 16.The proposed algorithm for MCCEKF-AKB (Figure 4) achieves a good estimation of the results with a fixed kernel bandwidth size.For instance, it is observed that among all the plots of EKF, MCCEKF, MCCEKF-AKB techniques, and the MCCEKF-AKB achieve the desirable performance and better accuracy with a fixed kernel bandwidth.Moreover, the suggested approach for determining the maximum kernel size is more reliable and effective.When compared to previous robust filters, the filter often maintains a better and more consistent estimation accuracy within the fixed range [14,15,20,27], even though the choice of confidence level variables affects the position of outliers.Consequently, MCCEKF-AKB's perfection and stability are verified.
good estimation of the results with a fixed kernel bandwidth size.For instance, it is observed that among all the plots of EKF, MCCEKF, MCCEKF-AKB techniques, and the MCCEKF-AKB achieve the desirable performance and better accuracy with a fixed kernel bandwidth.Moreover, the suggested approach for determining the maximum kernel size is more reliable and effective.When compared to previous robust filters, the filter often maintains a better and more consistent estimation accuracy within the fixed range [14,15,20,27], even though the choice of confidence level variables affects the position of outliers.Consequently, MCCEKF-AKB's perfection and stability are verified.bandwidth.Moreover, the suggested approach for determining the maximum kernel size is more reliable and effective.When compared to previous robust filters, the filter often maintains a better and more consistent estimation accuracy within the fixed range [14,15,20,27], even though the choice of confidence level variables affects the position of outliers.Consequently, MCCEKF-AKB's perfection and stability are verified.MCCEKF-AKB achieve the desirable performance and better accuracy with a fixed kernel bandwidth.Moreover, the suggested approach for determining the maximum kernel size is more reliable and effective.When compared to previous robust filters, the filter often maintains a better and more consistent estimation accuracy within the fixed range [14,15,20,27], even though the choice of confidence level variables affects the position of outliers.Consequently, MCCEKF-AKB's perfection and stability are verified.Finally, in Figure 17, the kernel size for the proposed study has been verified from the time history plot, which signifies the impact of kernel size.It can be seen that for both the Scenarios, MCCEKF-AKB shows qualitative performance in comparison to the other two methods, which also conveys the nature of the significant navigation preciseness efficiency.The aforementioned representations for kernel size in Figures 11 and 17 indicate the improved performance that MCCEKF updates the posterior estimates under MCC using a fixed-point iterative approach.The initial value of the iteration can be set to the previous estimate k x , and the tiny positive parameter offers a stop condition to acquire the convergent solution of the fixed-point iteration.It should be noted that when the kernel bandwidth, or  , in the MCCEKF is greater than a particular threshold, the conver- gence of the fixed-point algorithm under MCC is assured.Generally speaking, an algorithm's robustness and convergence speed increase with decreasing kernel bandwidth.However, the performance of the suggested algorithm will converge to that of the EKF algorithm as the kernel bandwidth increases, which verifies the consistent effectiveness of the navigation information processing [31].Finally, in Figure 17, the kernel size for the proposed study has been verified from the time history plot, which signifies the impact of kernel size.It can be seen that for both the Scenarios, MCCEKF-AKB shows qualitative performance in comparison to the other two methods, which also conveys the nature of the significant navigation preciseness efficiency.The aforementioned representations for kernel in Figures 11 and 17 indicate the improved performance that MCCEKF updates the posterior estimates under MCC using a fixed-point iterative approach.The initial value of the iteration can be set to the previous estimate x k , and the tiny positive parameter offers a stop condition to acquire the convergent solution of the fixed-point iteration.It should be noted that when the kernel bandwidth, or σ, in the MCCEKF is greater than a particular threshold, the convergence of the fixed-point algorithm under MCC is assured.Generally speaking, an algorithm's robustness and convergence speed increase with decreasing kernel bandwidth.However, the performance of the suggested algorithm will converge to that of the EKF algorithm as the kernel bandwidth increases, which verifies the consistent effectiveness of the navigation information processing [31].

Conclusions
In this paper, the algorithms for dealing with non-Gaussian or heavy-tailed errors, and outliers with multipath interference have been investigated by using three different filtering techniques.Since the multipath attenuation is crucial for improving positioning accuracy, a number of estimate techniques have been researched for high-precision positioning systems.
The MSE condition is used to create the standard EKF approach, which is constrained to the presumption that characterizes the linear and Gaussian ideal situation.Using an optimization strategy based on the correntropy criterion increases the resilience of nonlinear filters.A GPS navigation method based on MCC has been developed; it is constructive for interference errors with heavy tails and impulsive patterns.In addition, the behavior of the innovation-related variables has been implemented with the AKB for addi-

Conclusions
In this paper, the algorithms for dealing with non-Gaussian or heavy-tailed errors, and outliers with multipath interference have been investigated by using three different filtering techniques.Since the multipath attenuation is crucial for improving positioning accuracy, a number of estimate techniques have been researched for high-precision positioning systems.
The MSE condition is used to create the standard EKF approach, which is constrained to the presumption that characterizes the linear and Gaussian ideal situation.Using an optimization strategy based on the correntropy criterion increases the resilience of nonlinear filters.A GPS navigation method based on MCC has been developed; it is constructive for interference errors with heavy tails and impulsive patterns.In addition, the behavior of the innovation-related variables has been implemented with the AKB for additional performance enhancement, which is helpful in constructing the MCCEKF-AKB.
To demonstrate the effectiveness, simulations of GPS navigation tests have been performed for different filtering techniques.The adaptive filtering algorithm that employs the kernel bandwidth concept exhibits good potential as a substitute for the GPS navigation processor because it significantly improves navigation accuracy when compared to conventional methods, particularly the observations with non-Gaussian errors.For the purpose of illustration, two situations have been presented: (1) pseudorange observable errors based on Gaussian mixture distribution; and (2) pseudorange observable errors involved in the outlier type of multipath interferences with time-varying variance in measurement noise.The performance of other techniques, such as EKF, MCCEKF, and MCCEKF-AKB, have been compared, and the kernel bandwidth-based EKF algorithm, i.e., MCCEKF-AKB, has shown promising results in the enhancement of navigational accuracy.

Figure 2 .
Figure 2. The values of adaptive factor as a function of kernel bandwidth.

Figure 2 .
Figure 2. The values of adaptive factor as a function of kernel bandwidth.

Figure 2 .Figure 3 .
Figure 2. The values of adaptive factor as a function of kernel bandwidth.

Figure 3 .
Figure 3. Flow chart for calculation of the kernel bandwidth.

( 4 )
Computation of the measurement innovation based on MCC to obtain the factor: L k ; (5) Computation of the modified Kalman gain matrix: K k ; (6) Updating the state vector and the error covariance: xk and P k ; (7) Repeating from Step (2) to perform the subsequent estimation cycle.Sensors 2023, 23, x FOR PEER REVIEW 10 of 19

Figure 5 .
Figure 5. (a) The skyplot [19] and (b) the test trajectory for the simulated vehicle during the simulation interval.

Figure 5 .
Figure 5. (a) The skyplot [19] and (b) the test trajectory for the simulated vehicle during the simulation interval.

1  and 2 
represent the standard deviations of the individual Gaussian distri- bution, and the error model contamination perturbing parameter is  .To illustrate the viability of the suggested algorithm's robustness in the form of non-Gaussian distribution, the concerned requirements are 1 have been used to simulate the outlier type of multipath errors.The results of positioning errors for different techniques have been verified from Figures 7-10.The position errors corresponding to different values of extended kernel bandwidth factor  with respect to the time (s) are presented in Figure

Figure 6 .
Figure 6.Information for the sequence involving outliers interference for Scenario 1.

Figure 11 .
Figure 11.Time history of the kernel size for Scenario 1.

Figure 11 .
Figure 11.Time history of the kernel size for Scenario 1.

Figure 11 .
Figure 11.Time history of the kernel size for Scenario 1.

Figure 12 .
Figure 12.Information for the sequence involving outliers interference for Scenario 2.

Figure 12 .
Figure 12.Information for the sequence involving outliers interference for Scenario 2.

Figure 16 .
Figure 16.The results of position errors with three different techniques for Scenario 2 with β = 0.3: (a) East; (b) North; (c) Altitude.

Figure 17 .
Figure 17.Time history of the kernel size for Scenario 2.

Figure 17 .
Figure 17.Time history of the kernel size for Scenario 2.